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ABSTRACT 

The present study performs a six degree-of-freedom 
entry dispersion analysis for the Multiple Experiment Trans- 
porter to Earth Orbit and Return (METEOR) mission. 
METEOR offered the capability of flying a recoverable 
science package in a microgravity environment. However, 
since the Recovery Module has no active control system, an 
accurate determination of the splashdown position is diffi- 
cult because no opportunity exists to remove any errors. 
Hence, uncertainties in the initial conditions prior to deorbit 
burn initiation, during deorbit burn and exo-atmospheric 
coast phases, and during atmospheric flight impact the splash- 
down location. This investigation was undertaken to quanti- 
fy the impact of the various exo-atmospheric and atmo- 
spheric uncertainties. Additionally, a Monte-Carlo analysis 
was performed to statistically assess the splashdown disper- 
sion footprint caused by the multiple mission uncertainties. 
The Monte-Carlo analysis showed that a 3 -a splashdown 
dispersion footprint with axes of 43.3 nm (long), -33.5 nm 
(short), and 10.0 nm (crossrangc) can be constructed. A 58% 
probability exists that the Recovery Module will overshoot 
the nominal splashdown site. 

NOMENCLATURE 

a angle-of-attack (angle between the velocity vector and 
vehicle's xy plane), deg 

a T total angle between the velocity vector and the vehicle's 
axis of symmetry, deg 

P sideslip angle (angle between the velocity vector and 
vehicle’s xz plane), deg 
y flight-path angle, deg 
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INTRODUCTION 

The Multiple Experiment Transporter to Earth Orbit 
and Return (METEOR) mission was launched aboard the 
inaugural flight of the Conestoga launch vehicle in Octo- 
ber 1995. METEOR was a commercial program partially 
funded by NASA and developed by EER, Systems Incor- 
porated. It was formerly known as the COMmercial Ex- 
periment Transporter (COMET) with the goal of flying 
experiments to orbit in a microgravity environment and 



Fig. I Recovery module configuration. 


then be able to recover them. The spacecraft consisted of 
a Service Module (SM) which remains in orbit and a Re- 
covery Module (RM) which re-enters approximately 20 
days after launch. The RM was originally designed by 
Space Industries, Inc., and an entry dispersion analysis of 
the mission was performed. 1 *2,3,4 However, in 1995, EER 
took over the final development of the RM (Fig. 1), and 
modified the original mission (orbit altitude, spacecraft 
mass properties, and the nominal landing site). Also, the 
reliability of previous dispersion analysis was questioned 
due to a lack of credible aerodynamic data for the RM. 5 
As a result, a six degree-of-freedom (DOF) entry disper- 
sion analysis of the new mission scenario was necessitated. 

The mission scenario is as follows (Fig. 2): METEOR 
was to be launched into approximately a 250 nm circular 
orbit from NASA Wallops having an orbital inclination of 
40.5°. After approximately 20 days, the deorbit sequence 
begins. The attitude control system aboard the SM points the 
combined system to a specified attitude. With the use of a 
spin table, the RM is spun up to 73 rpm and separated from 
the SM. Since the RM has no active control system, the spin 
up provides stability to maintain its pointing attitude prior to 
performing the deorbit burn. Following a coast period (ap- 
proximately one half orbit), the deorbit burn is performed to 
initiate the entry of the RM. Shortly before the atmospheric 
passage (approximately 500 kft), the RM is despun to a nom- 
inal 6.5 rpm by a yo-yo despin device. This despin maneuver 
is performed to reduce the gyroscopic stability of the RM; 
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Fig. 2 METEOR mission profile. 


thus, allowing the vehicle’s aerodynamic characteristics to 
reduce its angle-of-attack prior to peak heating. The vehicle 
aerodynamically decelerates from approximately 25,000 ft/s 
to subsonic speeds. At an altitude of 60 kft, a series of three 
parachute deployments begin to slow the RM before splash- 
down off the coast of Virginia. The RM is positively buoy- 
ant and will be picked up by a surface vessel so that the ex- 
periments can be recovered. Unfortunately, due to a failure 
aboard the Conestoga during ascent, the launch vehicle and 
METEOR were lost. 

ANALYSIS 

Aerodynamics 

The aerodynamic database utilized in the flight simula- 
tion studies was derived from a combination of computa- 
tional fluid dynamics calculations, wind tunnel, and engi- 
neering code results. Since the flight traverses different flow 
regimes (rarefied, transitional, continuum), a range of solu- 
tion methods were required to estimate the aerodynamics of 
the RM. Free molecule models were employed at altitude 
above 394 kft (120 km); Direct Simulation Monte Carlo 
(DSMC) models 6 were used for altitudes between 295 kft 
(90 km) and 394 kft (120 km); and a Navier-Stokes solver 7 
was applied at altitudes below 263 kft (80 km). The Navicr- 
Stoke solutions were augmented with experimental data at 
Mach 6 and at Mach numbers between 1 .60 and 2.85. Sub- 
sonic aerodynamic characteristics were obtained with a lin- 
earized potential flow solver. In the transition region between 
263-295 kft (80-90 km), the aerodynamics are obtained 
through linear interpolation. The continuum aerodynamic 
characteristics of the METEOR Recovery Module are ex- 
plained in detail in Ref 8. Additionally, due to the RM\s 
configuration similarity with the Mercury manned entry cap- 
sule, dynamic stability aerodynamics were estimated based 
on Mercury-capsule flight data.^ 

Six DOF Trajectory' Simulation 

The trajectory analysis was performed using the six and 
three DOF versions of the Program to Optimize Simulated 
Trajectories (POST) 10 . This program has been utilized pre- 
viously in similar applications. 1 1 • ^ The six DOF version of 
this program was used to integrate the equations of motion 
from prior to the deorbil burn to parachute deployment. The 
three DOF program was then used from parachute deploy- 
ment to splashdown. The six DOF simulation includes the 
Earth atmospheric (GRAM-95) 13 and gravitational models, 
vehicle aerodynamics, mass properties and propulsion mod- 
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els, SM separation, and yo-yo release/despin models. The 
yo-yo release/despin is modeled using the method developed 
by Etter and Shamey. 1 ^ Their method models three-dimen- 
sional motion (i.e., precession and nutation), and allows for 
both center-of-gravity offsets and products of inertia. The 
three DOF simulation includes non-instantaneous parachute 
deployment models and the jettisoning of the RM’s aft end- 
cap. Since the RM has no active control system, no formal 
guidance and control strategy was utilized. 

Numerous sources of uncertainty affect the METEOR 
flight model. One source for this uncertainty arises through 
system capability limitations (e.g., the attitude control sys- 
tem on the SM can only achieve the prescribed pointing di- 
rection to within a specified tolerance). A lack of knowledge 
concerning the flight-day density, pressure, and winds and 
the computational uncertainty of the aerodynamics analysis 
are also contributing sources of error. Furthermore, measure- 
ment limitations in the mass, moments of inertia, center-of- 

Table i. Exo-Atmospheric Mission Uncertainties 

Mass Properties 

Mass 

eg position along spin axis .... 

eg position off spin axis 

Major moment of inertia 

(l X x- ’yy 

Cross products of inertia 

(l X y. Ixz- 'yz) 

Post-Separation State Vector 

Radial position 

In-track position 

Cross-track position 

Radial velocity 

In-track velocity 

Cross-track velocity 

Pitch/yaw attitude 

Pitch/yaw rate 

Roll rate 

Solid-Rocket Event 

Burn initiation time 

Rocket motor temperature .... 

Impulse 

Thrust vector cant angle 

Yo-Yo Event 

Despin initiation time 

Terminal roll rate 

Separation 

Spring induced velocity 


Table 2. Atmospheric Mission Uncertainties 


Atmospheric Flight 


T ransitional aerodynamics, C A 

+5% 

C N’ C Y 

+0.06 

Crm 

+0.01 

Continuum aerodynamics above Mach 10, 

C A 

±2% 

C N’ C Y 

±0.05 

Crrv C n 

+0.003 

Continuum aerodynamics below Mach 5, 

C A 

+10% 

C N' C Y 

+0.05 

^m 1 ^n 

±0.005 

Dynamic stability coefficients, C m q, C nr 

±50% 


Pressure, GRAM-95 model 3-o scale factors (3) 


Density, GRAM-95 model 3-o scale factors (3) 

Winds, GRAM-95 model 3-o scale factors (3, 

each component) 


Parachute deployment altitude (pilot/drogue) ±1000 ft 

Parachute deployment altitude (main) ±500 ft 

Parachute aerodynamics, C A ±20% 


gravity, and thrust nozzle alignment will also result in un- 
certainties. In this analysis, an attempt was made to conser- 
vatively quantify and model the degree of uncertainty in each 
mission parameter. 

For this mission, 57 potential uncertainties were identi- 
fied. These uncertainties are grouped into two categories (exo- 
atmosphcric and atmospheric) and are listed in Tables 1 and 
2, respectively, along with the corresponding 3-a variances. 
As seen, a few of the uncertainties have multiple entries (in 
parentheses) to account for variations at different mission 
phases. For example, in Table 1 , there are two center-ol-grav- 
ity off-axis position uncertainties (denoted by a 2 in the pa- 
rentheses). The first is applied prior to the deorbit burn to 
account for center-of-gravity measurement errors before 
launch, while the second is applied after the deorbit burn to 
account for unsymmctrical burning. For modeling the atmo- 
spheric properties (pressure, density, and winds: North-South, 
East-West, and vertical), three uncertainties are used to pro- 
vide a variation with altitude. Note the 3-a uncertainties list- 
ed in Tables l and 2 are estimates based on judgments of 
various experts in the field, and arc not based on statistical 
analysis (see Acknowledgments). 


±2 lb 

±0.25 in. (2) 

±0.25 in. (2) 

± 1 % 

±0.15 slug-ft^ 

±656 ft 
±8202 ft 
±1312 ft 
±0.66 ft/s 
±1.08 ft/s 
±1.12 tt/s 

} correlated 

±22 deg/sec 

±0.5 s 
±10°F 
±0.5% 

±0.45° 

±0.5 s 
±3.75 rpm 

0.226 ft/sec 
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RESULTS AND DISCUSSION 


Table 3. Nominal Mass Properties of the Recovery Module 


Nominal Mission 

Figures 3a and 3b show the groundtrack of the nominal 
entry trajectory, where the major mission events are highlight- 
ed. The splashdown point is at a latitude of 36.93° N and a 
longitude of 73.87° W, which is about 95 nm off the coast of 
Virginia. Table 3 lists the nominal set of mass properties of the 
RM. Note, in Table 3, the RM has a small nominal engine cant 
angle value; thus, the nominal engine thrust vector does not 
go through the center-of-gravity (eg) of the vehicle. 15 
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(b) Terminal mission segment. 
Fig. 3 Nominal mission groundtrack. 



Deorbit 

Entry 

Weight, lb 

763.0 

686.5 

Center-of gravity, ft: 



Along spin axis (x-direction, 



from nose) 

1.713 

1.597 

Off spin axis (y-direction) 

-0.00161 

-0.00161 

Off spin axis (z-direction) 

-0.00213 

-0.00213 

t X x» slug-ft 2 

46.40 

46.23 

lyy, Sllig-ft 2 

41.58 

40.43 

I zz , slug-ft 2 

40.47 

39.15 

I xy , slug-ft 2 

0.0191 

0.019 

l x2 , slug-ft 2 

-0.0191 

-0.0185 

I zz , slug-ft 2 

0.027 

0.0261 

Engine cant angle, deg 

0.22 

0.22 


Figures 4 and 5 show flight characteristics of the nomi- 
nal entry profile. The vehicle aerodynamically decelerates 
from approximately 25,000 ft/s to subsonic speed. During 
this time, the vehicle’s aerodynamic stability reduces the to- 
tal angle-of-attack to a moderate value (near 15°), as shown 
in Fig. 5. Peak heating occurs near Mach 2 1 with an angle of 
attack between 15° and 20°. The nominal attitude motion 
may be described by the superposition of two cycles, one in 
a, the other in (5, each of which is centered on 0°. Because 
the two cycles are out of phase, the total angle-of-attack, aj 
(defined as the angle between the velocity vector and the 
vehicle’s axis of symmetry), never reaches 0°. As shown in 
Figs. 4a and 5. a majority of the vehicle’s attitude decrease 
occurs while Hying in the transitional aerodynamic regime. 
This is the reason for the computational investment in the 
DSMC flowficld solutions. Without these computational so- 
lutions, the validity of this diminishing attitude motion (which 
is critical to mission success) would be questionable. 

Also shown in Fig. 5 is an increase in the total angle-of- 
attack towards the end of the entry. This behavior is an out- 
come of the sharp decrease in flight-path angle observed at 
the end of the entry (Fig. 4a). This area is also the region of 
Hight (lower velocities) in which dynamic-stability issues 
begin to dominate. If the vehicle was dynamically unstable, 
the RM would not be able to follow this flight-path angle drop 
and extremely high angles-of-attack would result. In the ex- 
treme case, a tumbling motion could result prior to parachute 
deployment. However, because the vehicle is dynamically 
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Fig. 4 Nominal mission profile. 
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Fig. 5 Nominal mission attitude profile. 


stable at all but the smallest total angles-of-aliack. a tumbling 
motion does not occur. Rather, the vehicle s dynamic-stabil- 
ity helps to minimize the effect of the llight-path angle vari- 
ations. Consequently, the total anglc-of-attack at parachute 
deployment for the nominal mission is on the order of 25°. 


Beginning at an altitude of 60 kft (approximately Mach 
0.8), a series of three parachutes are deployed which slow 
the RM down to approximately 19 ft/s prior to splashdown. 
Each of these parachute deployments is evident in Fig. 4b. 

Independent Uncertainty Effects 

To identify the mission uncertainties which have the 
greatest impact on the splashdown dispersion, each mission 
uncertainty was varied independently at its respective ±3-a 
value. Figure 6 shows the resulting splashdown range dis- 
persions for the largest contributors to the splashdown dis- 
persion size. Mission uncertainties not depicted in Fig. 6 had 
splashdown dispersions less than 2 nm. 


As can be seen from Fig. 6, the mission parameters have 
a varying effect on the splashdown dispersion. The disper- 
sions can be grouped into three categories: major (mission 
uncertainties 1-2), moderate (mission uncertainties 3-9), and 
small (mission uncertainties 10-20). The first group respon- 
sible for the largest dispersions, on the order of 25-30 nm, 
only include the exo-atmospheric uncertainties. These un- 
certainties produce the largest dispersion because they alter 
the deorbit burn direction. For example, a center-ol-gravity 
offset (mission uncertainty 1 ) causes a change in the loca- 
tion of where the thrust vector is applied from the nominal 
mission; producing the large splashdown dispersion. The 
second group of uncertainties produce dispersion on the or- 
der of 10-20 nm. Included in this group are several 
exo-atmospheric e fleets (specific impulse, rocket-nozzle cant- 
angle, and the initial weight uncertainties) and several 
atmospheric elYects t mid-altitude density and low-altitude 


Mitston Uncertainty 

1 Deorbii otl axis eg (0.25 m ) 

2 initial atirtude/rates 


Independent 

dispersion 

contribution. 

nm 


35 - 


30> 


25 ' 


20 - 

15- 


10 - 


5- 


| ♦ 7 n uncertainty 
Q -3 n unc«*n«mf> 



3 isp (0 5% t 

4 Rocxei nozzle cant angle (0 45°) 

5 initial weight (2 lb) 

£ M<d attitude density (GRAM 95) 

? Low altitude, east-wind (GRAM 95) 

8 Low attitude, north-wind (GRAM 95) 

9 Entry ofl axis C g (0 25 in ) 

i , Deorpit burn initiation (0 5 S) 
i i l K? (0 15 stug-tt 2 ) 

V Solid motor temperature (10°F) 

13 Atmospheric roll rate (5%) 
t4 Cross trac* velocity error (1.12 tVs) 

15 Radial velocity error (0 66 tt/s) 

16 Cross-trac* position error (1312 tt) 

17 I (0 15 siug-tl 2 ) 

18 Radial posmon error (656 ft) 

19 l KK MVi 

20 tyy (1%) 


man lemon 


2 3 4 5 t 7 8 9 10 11 12 13 14 1$ 16 17 18 19 20 
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Fig. 6 Significant contributors to the total splashdown 
range dispersion (3-0 variance shown in parenthesis).. 



horizontal and vertical winds uncertainties). The third group 
of uncertainties produce dispersions less than 5 nm. These 
uncertainties include those resulting from initial state-vector 
errors, moments of inertia misprediction, variation in the at- 
mospheric roll rate, and propulsion system unknowns 
(temperature and burn initiation time). 

As seen in Fig. 6, some of the mission uncertainties have 
a drastic difference in the dispersion between the +3-0 and 


Table 4. Major Contributors to Total Splashdown 
Range Dispersion 


Dispersion 
with +3 -<j 
U ncertainty 
(nm) 

Dispersion 
with -3-o 
Uncertainty 
(nm) 

1 . Deorbit off-axis eg (0.25 in.) 

15.5 

32.1 

2. Initial attitudes/rates 

29.0 

25.8 

3. Isp (0.5%) 

18.4 

17.6 

4. Rocket nozzle cant angle 

(0.45°) 

17.9 

12.2 

5. Initial weight (2 lb) 

9.7 

11.3 

6. Mid-altitude density 

(GRAM 95) 

9.9 

10.5 

7. Low-altitude, east-wind 

(GRAM 95) 

9.3 

10.2 

8. Low-altitude, north-wind 

(GRAM 95) 

9.3 

9.0 

9. Entry off-axis c.g. (0.25 in.) 

5.2 

9.2 

10. Deorbit burn initiation 

(0.5 sec) 

0.5 

4.9 

11. I xz (0.15 slug-ft 2 ) 

4.3 

3.7 

12. Solid-motor temperature 

(10°F) 

4.2 

0.6 

13. Atmospheric roil rate (5%).. 

2.3 

3.8 

14. Cross-track velocity error 

(1.12 ft/s) 

3.5 

3.7 

15. Radial velocity error 

(0.66 ft/s) 

0.3 

3.5 

16. Cross-track position error 

(1312 ft) 

2.6 

3.0 

17. I xy (0.15 slug-ft 2 ) 

0.8 

2.7 

18. Radial position error (656 ft) 

1 2.3 

2.6 

19. I xx ( 1% ) 

0.9 

2.5 

20. i yy (i%> 

2.4 

2.3 

Approximate 3-a Total 

53.6 

61.0 


All other contributors < 2 nm 


-3 -a values (mission uncertainties 1,9, 10, 12, 15, 17, 19). 
This outcome is due to the RM having a non-zero nominal 
engine cant angle value. Consequently, when the ±3-a val- 
ues are applied (for one of these mission uncertainties), one 
bound heightens the effect of the engine cant angle, while 
the other boundary value counteracts this effect. If there was 
no nominal engine cant angle, the resulting dispersions caused 
by the ±3-0 uncertainties would be more symmetric. 

The one-variable-at-a-time results are summarized in 
Table 4. Note that relative to the exo-atmospheric unknowns, 
the atmospheric uncertainties do not have a major impact on 
the splashdown dispersion. Furthermore, the aerodynamic 
uncertainty associated with the mission has a minimal im- 
pact on the splashdown dispersion size. Computing the L 2 
norm of these one-variable-at-a-time results, an upper-bound 
on the resulting 3-0 range dispersion from the Monte-Carlo 
analysis of no more than 50-60 nm is expected. 

Multiple Uncertainty Effects 

To determine the effects of multiple uncertainties occur- 
ring during the entry, a Monte-Carlo analysis is performed. 
Over 3500 random trajectories were simulated to assure a 
gaussian distribution. 

For some Monte-Carlo cases, a high amplitude oscillato- 
ry behavior in the total angle-of-attack near atmospheric in- 
terface was observed. As a result, for this vehicle, the total 
angle-of-attack at atmospheric interface is difficult to pin- 
point. This behavior can be observed in Fig. 7 (starting at an 
altitude of 500 kft ) for a particular Monte Carlo case, where 
the total angle-of-attack at atmospheric interface can be any- 
where between 60 - 1 20 . This phenomenon is a consequence 



Fig. 7 Oscillatory motion in total angle-of-attack. 
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of high pitch and yaw energies which are produced during 
the deorbit burn due to the thrust vector misalignment, cant 
angle, and the center-of-gravity uncertainty. During the ear- 
ly stages of the mission, these pitch and yaw energies are 
suppressed by the high roll rate (73 rpm) of the vehicle which 
provides inertial stability. With the RM despun to 6 rpm, roll 
stability is reduced and the pitch and yaw rates dominate the 
rotational motion. This behavior is not present in the nomi- 
nal mission (Fig. 5). However, a slight perturbation in one or 
a combination of the mission uncertainties can lead to this 
oscillatory motion. Approximately, 25% of the Monte-Carlo 
cases displayed this behavior. 

As a result of the oscillatory motion, the mean oscillato- 
ry value of the total angle-of-attack at atmospheric interface 
is used to express the true attitude of the RM. Figures 8-10 
show the distribution of the total angle-of-attack at atmo- 
spheric interface, peak heating, and parachute deployment. 
At atmospheric interface, the mean of the mean total angle- 
of-attack is 66.8° with an oscillation amplitude of 17.2°. At 
peak heating and parachute deployment, the mean total an- 
gle-of-attack is 20.6° and 5 1 .2°, respectively. The maximum 
total angle-of-attack at each of these events can be signifi- 
cantly higher than these mean values. Maximum values of 
1 76.6° ( 1 1 1 .4° plus a 65.2° oscillation) at atmospheric inter- 
face, 56.7° at peak heating, and 170.5° at parachute deploy- 
ment were obtained. Table 5 summarizes these results. Note, 
the heatshield design limit on the total angle-of-attack at at- 
mospheric interface of 75° is v iolated in some cases. How- 
ever. since the frequency of the oscillatory motion is high, 
the heating levels on the side and aft end of the RM (for 
these high total angle-ol-altack cases) were determined to he 
acceptable. 16 * 17 In addition, the reliability of the aerodynamic 

1200 r- , 

! HI' ! 


1000 - 



a*p. deg 


Fig. 8 Distribution of the mean total angle-of-attack at 
atmospheric interface resulting from over 3500 
Monte-Carlo simulation cases. 



Fig . 9 Distribution of the total angle-of-attack at peak 
heating resulting from over 3500 Monte-Carlo 
simulation cases . 



uj. deg 

Fig. 10 Distribution of the total angle-of-attack at 
parachute deployment resulting from over 3500 
Monte-Carlo simulation cases. 

database for high altitudes and angle-of-attaeks greater than 
145° is questionable, because these conditions exceed the 
limit of the database. Approximately. 0.39} of the Monte- 
Carlo cases had total angle-of-attaeks greater than 145.0 

Figures I 1 and 12 show the downrange and crossrange 
distribution at splashdown, respectively. The minimum down- 
range is -43.8 nm (short) from the nominal splashdown point, 
while the maximum downrange is 67.0 nm (long). The max- 
imum crossrangc obtained is I 1 .5 nm. A 3-0 dispersion foot- 
print with axes of 43.3 nm (long). -33.5 nm (short), and 10.0 
nm (crossrangc) can be constructed. Within the assumptions 
of the present analysis, a 99.739} probability exists that the 
RM will splashdown within this 3-0 footprint. Table 5 sum- 
marizes these and other statistical results. 
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Table 5. Summary of Monte- Carlo Analysis 



Mean 

Min 

Max 

3-c 

Mean atmospheric interface attitude, deg ... 

.... 66.8 

36.6 

111.4 

26.2 

Amplitude about mean atmospheric 
interface attitude, deg 

.... 17.2 

1.6 

65.2 

28.7 

Peak heating attitude, deg 

.... 20.6 

8.5 

56.7 

12.8 

Parachute deployment attitude, deg 

51.4 

14.9 

170.5 

56.7 

Splashdown downrange, nm 

2.6 

-43.8 

67.0 

43.3 (long) 

Splashdown crossrange, nm 

-0.5 

-11.5 

10.9 

-33.5 (short) 
10.0 

Total splashdown range, nm 

11.0 

0.2 

67.6 

23.8 

Splashdown latitude, deg 

... 36.83 

36.52 

37.17 

0.26 

Splashdown longitude, deg 

... 286.10 

285.22 

287.40 

0.76 



aj. deg 

Fig. II Downrange at splashdown. 
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Fig. 12 Cross range at splashdown. 


Figures 13 and 14 show the atmospheric interface and 
splashdown locations for each of the cases. Note, the total 
downrange dispersion shown in these figures is similar. This 
outcome is due to the dominant effect which the exo-atmo- 
spheric uncertainties have upon downrange. While impact- 
ing downrange to some extent, the atmospheric uncertainties 
have a more pronounced effect on crossrange. Additionally, 
there is a slightly higher probability (58%) that the trajectory 
will splashdown long than short. 

The present results predict a larger splashdown disper- 
sion than previous METEOR entry dispersion analyses that 
were performed early in the program.^* 4 Table 6 summariz- 
es the findings from Refs. 3 and 4. The splashdown predic- 
tions from the present analysis arc approximately double that 
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Fig. 13 Range dispersion at atmospheric interface 
resulting from over 3500 Monte-Carlo simulation cases. 



of Refs. 3 and 4. The differences in these results can be at- 
tributed to many factors. The initial orbit altitude, target land- 
ing site, and mass properties of the RM have changed. More 
mission uncertainties are considered in the present analysis 
than in the previous studies. These additional uncertainties 
add to the overall dispersion. Additionally, the 3-0 variances 
selected in the present study are based on more conservative 
estimates (roughly twice that of the previous studies in some 
cases, e.g., Isp, solid-motor temperature, weight, and initial 
attitude/rates). Moreover, the aerodynamic characteristics of 
the RM are much better known due to the extensive compu- 
tational and experimental investment. Overall, the present 
analysis produced a more conservative estimate of the splash- 
down dispersion. 

SUMMARY 

The present study performs a six degree-of-freedom en- 
try dispersion analysis for the Multiple Experiment Trans- 
porter to Earth Orbit and Return (METEOR) Recovery Mod- 
ule. For this mission, 57 potential exo-atmospheric and 
atmospheric uncertainties were identified. From a one-vari- 
able-at-a-time uncertainty analysis, where each variable was 
set at its estimated ±3-0 value, a center-of-gravity offset from 
the spin axis and initial attitude/rate uncertainties where 
shown to produce the greatest dispersions (each on the order 
of 30 nm). Uncertainties in specific impulse, engine cant 
angle, initial weight, mid-altitude density, and low-altitude 
winds produced dispersions on the order of 10-20 nm each. 
All other uncertainties produced dispersions less than 5 nm. 
From a Monte-Carlo analysis of over 3500 random, off-nom- 
inal trajectories, a 3-0 splashdown dispersion footprint with 
axes of 43. 3 nm (long), -33.5 nm (short), and 10.0 nm (cross- 
range) was obtained. Within the assumptions of the present 
analysis, a 99.73W probability exists that the RM will splash- 
down within this 3-0 footprint. Additionally, there is a 589? 
probability that the Recovery Module will overshoot the nom- 
inal splashdown site. Furthermore, the present analysis pre- 
dict a larger splushdow n dispersion than previous METEOR 
entry dispersion analyses. 


Monte-Carlo Results 

22.4 (long) 17.0 (long) 
22.4 (short) -17.0 (short) 
3.8 2.9 


Downrange, nm 
Crossrange, nm 



Fig. 14 Range dispersion at splashdown resulting from 
over 3500 Monte-Carlo simulation cases . 


Table 6. Summary of Previous Entry Dispersion Analysis 


One-At-A-Time Results 



Space 

Aerospace 

Mission 

Industry 3 

Corp . 4 

Uncertainty 

Dispersion 

Dispersion 


(nm) 

(nm) 

Atmosphere (20% 3 , 1 0% 4 ) 

14.0 

8.5 

Isp (0.25% 3 - 4 ) 

10.3 

13.5 

Ballistic coefficient (5.0 Ib/ft^ 4 ) 

NA 

9.6 

Drag coefficient (10% 3 ) 

6.2 

NA 

Atmospheric roll rate (3.75 rpm 3 ) . 

6.3 

0 

Initial attitude (0.7 deg 3 , 0.5 deg 4 ) 

6.2 

0.8 

Initial weight (1.2 lb 3 ) 

6.0 

0 

State vector 

5.9 

3.2 

Winds 

Deorbit burn initiation 

4.2 

0 

(0.65 s 3 , 1 .0 s 4 ) 

2.4 

0.2 

Parachutes 

2.0 

0 

Solid-motor temperature (5°F 3 ) .... 

2.0 

0 

Thrust cosine loss (1.5 deg 3 ) 

1.0 

0 
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